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ABSTRACT 

This  note  develops  the  basic  B-spline  theory  without  using  divided  differences.  In¬ 
stead,  the  starting  point  is  the  definition  of  B-splines  via  recurrence  relations.  This  ap¬ 
proach  yields  very  simple  derivations  of  basic  properties  of  spline  functions  and  algorithms.^ 
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SIGNIFICANCE  AND  EXPLANATION 

The  standard  definition  of  the  univariate  B-spline  involves  the  divided  difference  and 
the  calculus  of  divided  differences  has  been  heavily  used  in  developing  the  univariate  spline 
theory. 

In  this  report  we  discuss  an  alternative  approach.  We  choose  the  recurrence  relations 
for  B-splines  as  the  basic  concept.  This  leads  to  very  simple  derivations  of  the  main 
properties  of  spline  functions  as  well  as  of  the  basic  algorithms. 
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B-SPL1NES  WITHOUT  DIVIDED  DIFFERENCES 


Carl  de  Boor11’  and  Klaus  Hollig  1  * 


1.  Basic  Properties 


Let  t  :  ••  •.  f_i.  be  a  nondecreasing,  biinfinite  sequence  of  real  numbers 

with  lirn,  .  -  x/,  ~  x.  The  B-splines  corresponding  to  the  “knot  sequence”  t  are 

defined  by  the  recurrence  relation 


Bi.k  ■  -  ^'i.k^i.k  1  (1  ^t-r  l.k)Bt .» J,fc-  1 


with 

BtA*) 

«.*(/ ) 

This  gives  B, >  in  the  form 


if  t,  t  <  t, . , 
ot  her  wise 

,  •  if  t,  <  C-fc-  l 
otherwise 


(la) 


(16) 


i  *  k  1 

X  (2) 

j 

with  each  bjk  a  polynomial  of  degree  k  since  it  is  the  sum  of  products  of  k  1  linear 
polynomials.  From  this  we  read  off  that  Bx_k  is  a  piecewise  polynomial  of  degree  <  k 
which  vanishes  outside  the  interval  t,.l,.k  and  has  possible  breakpoints  /,.•••, In 
particular.  is  just  the  zero  function  in  case  /,  Also,  by  induction.  Bt,k  is 

positive  on  the  open  interval  (f,./,  *  *).  since  both  a.',.*  and  1  are  positive  there. 

But  it  is  a  bit  of  a  miracle  that  the  recurrence  relations  produce  smooth  functions.  We 
discuss  this  question  in  Theorem  2  below. 

A  spline  of  order  k  with  knot  sequence  t  is.  by  definition,  a  linear  combination 
of  the  B-spiines  B,  k  associated  with  that  knot  sequence.  Let 

-S.t  :  {  X  a’B,.k  ■  <>,  R }  («0 

T 

denote  the  collection  of  all  such  splines.  We  now  explore  this  space. 

(1)  Sponsored  by  the  United  States  Army  unde  r  Contract  No.  DA A(i29-X0-C-(X)-1 1 . 

(2)  >upported  by  International  Business  Machines  Corporation  and  National  Science 
Foundation  Grant  No.  DMS-S351 1*7. 


We  deduct*  from  the  recurrence  relation  that 


\  •  «,  i ( 1 


•*'t.k))Bi.k  1  • 


On  t  tie  other  hand,  arguing  a*-  in  2  .  for  the  special  sequence 


(4) 


",  :  v,.fc(')  :  (0-t  k  i  r) 

(with  r  R).  we  find  for  /jf , . ^  i  -  0.  i.e..  for  /,  •  0-a  i  'hat 

o,  i(l  *,i)  -  i(r)  ((0 - 1  -  r)~Y*  *  (0  -  r)U  - 

t-'i.tc  l  (r)(-  ?} 

since  f  {t t  -  k  i )  — ' j . a-  ■  /(/,)( 1  is  the  straight  line  which  agrees  with  /  at  <,-*  *-  t  and 

t , .  This  shows  that 


y"f  v, .k{r)tti.k  -  (•  T)  Vt.k  j(r)B,,*-  j. 

hence,  by  induction,  that 

y^v,.k(')Bt.k  --  (•  -  r)k  1  ^V«.t(r)B,.,. 

This  proves  the  following  identity  due  to  Vlarsden. 

Theorem  1.  For  any  t  ■  R. 

(•  r)k  1  ,  y^r,,k(r)Bt.k.  (5) 

t 


with  r,.A-(r)  :  (0  .  ,  ')•••  (0  -At-  *)• 

Since  r  in  ("»|  is  arbitrary  it  follows  that  .S'a  t  contains  all  polynomials  of  degree  <  k. 
More  than  that,  we  can  even  give  an  explicit  expression  for  the  required  coefficients,  as 
follows. 

My  differentiating  (■">)  with  respect  to  r.  we  obtain  the  identities 


l>) 


k  I  i 


l'i.k(<  I 


{k  I)' 


If ,  A  •  1/  . 


(6) 


with  Df  the  derivative  o|  the  function  /  On  using  this  identity  in  the  Taylor  formula 


for  a  polynomial  p  of  degree  <  k ,  we  conclude  that  any  such  polynomial  can  be  written  in 
the  form 

p  =■■  *>.kP  Ri.k-  (") 


with  A,  >  given  by  the  rule 


w--  x; 


v-  LP^y.Mr, 

h,  (*■-')! 


Here  are  two  special  cases  of  particular  interest.  For  p  -  ].  we  get 

i  =  £>,.*  (9) 

since  Dk~lu\.k  -  (~l)k~l(k  -  1)!,  and  this  show’s  that  the  B,  *  form  a  partition  of 
unity.  Further,  since  Dk~2v t.k  is  a  linear  polynomial  which  vanishes  at  tj  (/,*  ,-•••- 
i)  (*  -  1).  _ 

f  -  ^f( t[)Bt<k  for  every  linear  polynomial  L.  (10) 

1 

The  identity  (5)  also  gives  us  various  piecewise  polynomials  contained  in  t:  If 
r  -  tj,  then  v,,k(T)  vanishes  for  i  -  j  -  k  ~  1.  •  ■  • ,  j  -  1.  Since  B,.*(/)  =  0  for  i  <  j  -  k  ~  ] 
and  t  <  tv  it  follows  that 

CX 

(-  -  (ii) 

«=j 

with  o. :=  max{o,0}  the  positive  part  of  the  number  q.  The  same  observation  applied 
to  (6)  shows  that 

(•-?;)*  €  Sk,  t  for  I  <  H  <  Hlj  *  {<,:#,=  *,}.  (12) 

Theorem  2.  If  f,  <  li  +  fc'  then  the  B-splines  B,  ^  are  linearly  independent  and 
the  space  Sk.t  coincides  with  the  space  S  of  all  piecewise  polynomials  of  degree  <  k  with 
breakpoints  t,  which  are  k  -  1  -  #f,  times  continuously  differentiable  at 

Proof.  It  is  sufficient  to  prove  that,  for  any  finite  interval  /  a.b  .  the  restriction 
S  j  of  the  space  S  to  the  interval  /  coincides  with  the  restriction  of  Sk, t  to  that  interval. 
The  latter  space  is  spanned  by  all  the  B-splines  having  some  support  in  /.  i.e..  all  B,  t 
with  (t,.t,.k)  /  ~  0.  The  space  5i/  has  a  basis  consisting  of  the  functions 

( •  a ) k  1  .  i>  1 .  •  •  ■ ,  k\  (•  -  t , )  * ~  u .  1/  -  1 .  •  •  • .  a  t , ,  for  a  <-  / ,  •  b.  ( 1  ) 

Tfiis  follows  from  the  observation  that  a  piecewise  polynomial  function  /  with  a  breakpoint 
at  t.  which  is  k  1  at,  times  continuously  differentiable  there  can  be  written  uniquely 

a  t 

J  ~  p  -  5l°'  (-  '■)*  "• 


as 


with  p  a  suitable  polynomial  of  degree  •  k  and  suitable  coefficients  a,  .  Since  each  of  the 
functions  in  (13)  lies  in  S'*.,,  by  (6)  and  (12).  we  conclude  that 

S  /  (SYt )  /•  (H) 

On  the  other  hand,  the  dimension  of  .S  ;.  i.e..  the  number  of  functions  in  (13).  equals  the 

number  of  B-spiines  with  some  support  in  /  (since  it  equals  k  -  ^2a<t  ,  rOi).  hence  is  an 

upper  bound  on  the  dimension  of  (S*  t)  /•  This  implies  that  equality  must  hold  in  (14), 
and  that  the  set  of  B-spiines  having  some  support  in  7  must  be  linearly  independent  over 
7. 

Corollary  1.  All  B-splines  having  some  support  on  a  given  interval  are  linearly 
independent  over  that  interval. 

Corollary  2.  If  t  is  a  refinement  of  the  knot  sequence  t.  then  5*ft  C  Sk  j. 

Corollary  3.  If  t,  <  (,**_).  then  the  derivative  of  a  spline  in  S*. t  is  a  spline  of 

degree  ■-  k  1  with  respect  to  the  same  knot  sequence,  i.e..  DSk,t  - 

The  identity  (7)  can  be  extended  to  all  spline  functions.  For  this,  we  agree,  consistent 
with  (lb),  that  all  derivatives  in  (8)  are  to  be  taken  as  limits  from  the  right  in  case  r 
coincides  with  a  knot. 

Theorem  3.  If  r  in  definition  (8)  of  A,  *  is  chosen  in  the  interval  /,.*,  +  *).  then 

A,.*( jBj.it)  =  Ox-  (I®) 

i 

It  is  remarkable  that  r  can  be  chosen  arbitrarily  in  the  interval  f,./,,*)-  The  reason 
behind  this  is  that  A,.*/  does  not  depend  on  t  at  all  if  /  is  a  polynomial  of  degree  less 
t  hail  k. 

Proof.  Assume  that  r  t)  :  /,.<,-*■)  and  let  p;  be  the  polynomial  which 

agrees  with  77,  *  on  (li.ti.  i).  Then 

^i.kBj.k  -  A x  ^Pj. 

On  the-  other  hand. 

i 

Pj  ^i.kP)  Pi- 

i  —  /—  I  k 

sine e  this  holds  by  (7)  on  f  j.fj .  i ).  while,  by  Corollary  1  or  directly  from  (7 ).  p; .  i  *.•••.  p/ 
are  linearly  independent.  Therefore  necessarily  X,.kPj  equals  1  if  i  /  and  0  otherwise. 


2.  Algorithms 

In  this  section,  the  basic  algorithms  for  computing  with  the  B-spline  representation 
are  derived.  Unless  otherwise  staled,  all  algorithms  refer  to  the  spline  function 

.s  y  <‘,1*1. k-  i  H») 

i 

In  the  following,  the  subscript  k  is  omitted  whenever  possible,  e.g.  B,  :  -  B,k-  v,  :  \',.k- 

etc. 

The  spline  s  can  be  evaluated  using  the  recurrence  relation.  By  (4). 

y  ° i f^i.k  —  y  ^  .k^ t  (i ~  ■‘m .k i - 1  j  b, ^  i  —  •  y  ° ,  b,  k _  i . 


Iterating  this  identity  one  finally  arrives  at 

*  -=  XX  'B'' 

1 

and.  by  the  definition  of  B, j,  the  right  hand  side  equals  a*“  1  on  t}.t j).  This  yields 

Algorithm  1  .  On  the  interval  s  -  a with  the  polynomials  a\  com¬ 

puted  as  follows: 


<*’  :  -  (1  -  ,K_ J  -  k  -  r  -  1  <  i  <.  j. 

For  the  relevant  range  of  indices.  ~'uk  r(0  -  ;0.  1  so  that  s(f)  is  computed  by  repeat¬ 
edly  forming  convex  combinations  of  the  B-spline  coefficients. 

It  follows  from  Theorem  2  that  the  derivative  of  .s  is  a  spline  of  degree  *  k  1  with 
respect  to  the  same  knot  sequence,  i.e. 

Ds  ,.  ( 1 7 1 

1 

By  Theorem  3, 

a',  -  A,.*.  ,(/As)  (18) 

if  r  is  chosen  in  the  interval  {t ,.  t,  _  k- i  )•  To  relate  a'  to  «,  we  express  A,>  , 1)  as  a  linear 
combination  of  the  functionals  X,  k.  making  use  of  the  fact  that  X,k  depends  linearly 
on  v,.k  and  that 

U<-k-  i  i  t’..t  <•'.  i  t-  (19) 

From  the  definition  (8). 


A,,  ,a/M  -  v»  (jt  ;■ 


\*i  *»*  a  ‘  tv*  1  ( r ) I}‘  * 1  /(7 )  -  v( 


-■fl'VW. 


the  last  oqualit\  b\  setting  p  :  ■  1  and  using  the  fact  that  Dr  t  0-  Comparison 

of  these  two  lines  shows  with  the  aid  of  (19)  that 

A<  *  'D  r  t  (x'  k  A'-1  *)• 

Assuming  that  Bt,k  t  *  0.  i.e..  that  /,  •  we  can  choose  r  t  = 

(/,  i-t,-k  i)  {tfti-k)-  B>  Theorem  3.  this  yields 

Algorithm  2.  Compute  the  coefficients  for  ^2a',Bi.k  -  \  D\^a,Bx_k  by 
'  a,  a,  i 

a  —  V V 777  <  M-t-fc-l* 

(/« .  *  1  -  /,)  (k  --  1) 

By  Corollary  3.  St  ;  5j  for  any  refinement  t  of  the  knot  sequence  t.  and  therefore 
any  spline  s  •  can  be  written  as  a  linear  combination  '£2alB,  of  the  B-splines  Bt  which 
correspond  to  the  refined  knot  sequence.  The  computation  of  the  new  coefficients  a,  from 
the  a,  constitutes  the  knot  insertion  or  subdivision  algorithm  used  in  CAGD  1.4  . 
For  this,  we  need  to  express  d,  in  terms  of  the  a,.  Bv  Theorem  3.  this  is  equivalent  to 
comparing  the  corresponding  A,  with  At.  Since  A,  depends  linearly  on  17, ,  this  requires 
nothing  more  than  to  express 

Vt  ~  (G  +  i  ~  ■)  •  •  •  (t  t-1- *  - 1  ■■  ■) 
a1-  a  linear  combination  of  the  w 

This  is  particularly  easy  when  t  is  obtained  from  t  by  adding  just  one  knot,  say  the 
point  t  t , .  1  j .  | )  so  that 

|  <>•  if  i  <  J- 
t ,  =-  |  t.  if  i  ~  j  ~  1 ; 

\  t,  i.  if  i  >  j  ->■  1 . 

Then 

'  J  k~  2; 

!  \  V,-  1-  '  J- 

hence  there  is  some  actual  computing  necessary  only  for  j  k  -  2  <  i  <  j.  For  this  case, 
observe  t  hat  / , .  (  •  t  I , , )  •  t , .  *  i  •  hence 

of,  •  ,h’t  [  (/, .  i  ■)■■■(/,.  k  "  •)  «(/,- k  i  ■  •)  *  ■!(>>  •) 


foyi 

’-1 

•  •  •  %  -  4 
•*.  *'.  •  1 


c.y.y.j 

yCVvO 


■-S  ■ 


‘‘  *v  .*• 


-6- 


provided  <>(<,-*-1  -  •)  -  3(t ,  -•)  =  (/-  •).  i.e.. 

a  =  «c,(/)  and  J  (1  -  -,(/)). 

Since/",.]  ■  /  •  wo  can  choose  r  in  the  definition  (8)  in  the  interval  (/,,/, 

|/,  ,)  (t,.t,-k).  This  proves 

Algorithm  3.  If  the  knot  sequence  t  is  obtained  from  the  knot  sequence  t  by 
addition  of  the  point  (  -l  tj.tj ,j).  then  the  coefficients  a,  for  the  spline  s  with  respect  to 
the  refined  knot  sequence  are  given  by 

if  1  ■  j  k  -  2: 

(1  -*,(£))<!, „1.  if  j  ~  k  ~  2  <  t  <  j:  (20) 

if  1  >  J- 

Observe  that  u.’(f)  c  [0, 1  for  the  relevant  range  of  indices  and  thus  the  coefficients  a 
are  convex  combinations  of  the  coefficients  a. 

If  r  <  Ar  -  1,  then,  after  just  (k  -  1  -  r)-fold  insertion  of  i,  we  obtain  a  knot 

sequence  t  in  which  the  number  t  occurs  exactly  k  -  1  times.  This  means  that  there  is 
exactly  one  B-spline  for  that  knot  sequence  which  is  not  zero  at  t.  Hence  it  must  equal  1 
at  i  and  its  coefficient  must  provide  the  value  of  s  at  t.  This  makes  it  less  surprising  that 
the  calculations  in  Algorithms  1  and  3  are  identical. 


|  a,. 

a,  -  ,  ^.’,(f)  a, 

I  a,  ]. 


V 


yjgr 


!ir 
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